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The first and second famihes of D-wave bb quarkonium states are expected 
to have masses near 10.16 and 10.44 GeV/c^. The accuracy of these 
predictions is discussed, and the prospects of two methods for produc- 
ing these states in electron-positron collisions are updated. Direct scans 
in the e+e~ center-of-mass can give rise to the ^Di states. The l^Dj 
states have also been searched for in the electromagnetic cascades from 
T{3S) •yXb ~^ 11^ D J. The sample of T(3S') decays required to defini- 
tively observe the l^Dj states is found to be only somewhat greater than 
the present world's total. 



PACS Categories: 14.40. Gx, 13.20.Gd, 13.40.Hq, 12.39.Ki 

The bound states of heavy quarks have provided a key testing ground for quantum 
chromodynamics (QCD). Perturbative QCD describes the short-distance aspects of 
the interquark force and many of the decays of the states, while non-perturbative 
aspects are probed by the long-distance interaction and by details of fine- structure and 
hyperfine splittings 

The ^5*1 states of charmonium (cc) and bottomonium {bb) are easily accessible in 
both hadronic Drell-Yan processes and electron-positron collisions, since they couple 
directly to virtual photons. These states then decay electromagnetically to others 
with sizeable branching ratios. In this way both the lowest charmonium state, the 
Tjc = cc(^S'o), and many P-wave cc and bb states have been discovered. With knowledge 
of the masses of the lowest Xc = cci^Pj) and the two lowest Xb = bb{^Pj) families, the 
interquark interaction can be mapped out to an extent which permits anticipation of 
other, as yet unseen, levels. Foremost among these levels are the rj^ = bb(^So) states, for 
which we have recently suggested observation strategies 0, and the D-wave bb levels. 
A candidate for the lowest ^Di cc level is V'lSTTO), which couples sufficiently strongly 
to e+e~ to be useful as a copious source of charmed meson pairs. 

In this note we review some predictions of bb{n^Dj) masses for n = 1 and n = 2, 
where n is the level number. These levels are expected to lie below the BI3 flavor 
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Figure 1: Predictions for the spin-weighted averages of ^Dj bb states. Open circle 
(o): KR 1^ (inverse scattering); open square (□): Cornell (QCD-based potential); 
open triangle (A): BT (QCD-based potential; published masses only quoted to 
nearest 10 MeV); open inverted triangle (y): Gupta et al. (QCD-based potential); 
open diamond (O): MR (QCD-based potential); solid circle: MB |T^ (relativistic 
corrections; mass oi^Di plotted); solid square: GI (QCD-based potential, masses 
calculated to nearest MeV); solid triangle: Grant et al. (power-law potential); 
solid inverted triangle: EQ (QCD-based potential [0, quoted for n = 1): solid 
diamond: lattice [T^ (quenched approximation with (3 = 6.0, quoted for n = 1). 



threshold and thus to be quite narrow. A comprehensive treatment of them was pre- 
sented in Ref. 0. We compare the predictions of that work with others, estimate the 
likely errors in predictions of masses, and update the prospects for discovering some 
of the predicted levels. These questions have taken on renewed interest as a result of 
plans by the CLEO Collaboration 0] to examine some aspects of T spectroscopy, both 
via significant augmentation of the world's sample of T{3S) decays and via direct scans 
for ^Di states in the 10.13-10.17 and 10.42-10.46 GeV/c^ mass range. 

We shall show that most potential models predict a narrow range of values for the D- 
wave masses. The question has been raised of whether any potential description is valid 
for quarkonium [Q; discovery of the D-wave states within the predicted range would 
not lay such doubts to rest, but would be one further point in favor of a description 
which, at least for the bb system, has been remarkably successful. 

Some predictions of spin-weighted masses of the n^Dj bb levels 0, ^ 0, ||, |^, |10 



11, 12, 13, [T^ are summarized in Fig. |T| [^. Most potentials give a center-of-gravity 
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Table I: Deviations of predicted masses for bb levels from observed values, in MeV/c^. 



Reference T(15) XbAlP) T(25) XbA2P) T(35) T(45) 



i (a) 3 

i (a) 23 

(b) (a) -10 

1 2 
i (a) 6 
(a) 23 
O (d) 5 -16 
(e) -5 -9 
4 -17 
113 (a) -13 ±17 
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40 


-20 


-8 


-1 


55 


19 


-2 
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-7 


-16 


-29 


-16 


22 
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105 ± 40 


30 ±80 


(c) 


10023 


10260 


10355 


10580 



(a) Input, (b) Published masses quoted to nearest 10 MeV. (c) Not quoted, 
(d) Numbers in this table are based on masses calculated to 1 MeV, while 
those published in Ref. 
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were rounded to the nearest 10 MeV. 

(e) Ref. quotes spin-averaged nS masses; values here are for states 

(f) Experimental masses quoted to nearest MeV/c^. 



of the ID levels within ±10 MeV/c^ of 10160 MeV/c^, except for the more recent 
treatment of Eichten and Quigg , based on the Buchmiiller-Tye [0 potential, which 
gives a value about 33 MeV/c^ lower. However, their model also underestimates several 
other bb masses in comparison with experimental values |T^, as shown in Table |. A 
quenched lattice calculation |l^ gives a l^D2-l^Si splitting of 761(20) MeV/c^ for 
P = 6/gQCD = 6.0 if the IP-IS splitting is used to set the scale. This becomes 720(30) 
MeV/c^ if the 2S'-1S' splitting is used to set the scale (see [0 in which it is pointed 
out that this scale is the one that makes the bare quark mass used closest to the b). 
The result with the 2S'-1S' scale is quoted in Fig. |l], resulting in a ^Dj mass of 10180 
MeV if one assumes a small singlet-triplet splitting in the ID levels. 

The difference of scales obtained using different experimental splittings is a fea- 
ture of quenched lattice simulations. This difference should disappear once physical 
dynamical calculations can be done. In the meantime it is an additional source of 
uncertainty. The most recent unquenched calculations are given in Ref. but clear 
signs of the effect of dynamical quarks are hard to see. For S-wave and P-wave lattice 
predictions, published results in the quenched approximation are given in Ref. [p^j . 
These are quoted in Tables | and |T| using the f3 = 6.0 results and setting the scale from 
the 25-15' splitting, for consistency with the results quoted in Fig. [l| |0. 

The predicted centers-of-gravity of the 2D levels range from 10430 to 10455 MeV/c^. 
These levels are more likely to be affected by coupling to B^*^B^*^ systems above 
the flavor threshold of 10560 MeV/c^. The fine-structure splittings predicted in various 
models are summarized in Table 0. 

The production of the lowest D-wave bb states is most likely in cascade transitions 
from the T(3S'). In Ref. [0], two means of studying these transitions were proposed. 
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Table II: Fine-structure splittings for spin-triplet states predicted in various models. 
Shown are deviations from spin- weighted centers-of-gravity, in MeV/c^. 



Reference 
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2D 







L - 1 


-40.4 (a) 
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(a) Input, based on experimental masses, (b) Published masses quoted to nearest 10 MeV; 
numbers in this table calculated to 1 MeV. (c) Not quoted. 

One can employ the three-photon inclusive transitions 35 ^ 2P ^ ID ^ IP, or 
the four-photon transitions 35* 2P ID IP — * 15* — > which should 

have considerably less background. Either process suffers from backgrounds due to 
3S ^2P ^2S ^ IP and 35 ^ lSn%^. In four-photon processes one can eliminate 
events in which two pairs of photons are consistent in mass with two neutral pions. 
However, in such processes the cascade 3S -+ 2P -+ 2S ^ IP ^ IS ^ l+l- leads 
to events in which the photons from 2P ^ 25" are easily confused with those from 
ID ^ IP, while those from 2S ^ IP are easily confused with those from 2P ^ ID. 

Using the branching ratios predicted in Ref. for the electromagnetic transitions 
from one bb state to another, and the measured branching ratio B{T{1S) —> e~^e~) = 
(2.38 ± 0.11)% [|^, one predicts the numbers of 47e"'"e~ events per 10^ T{3S) shown 
in Table |T|. The total numbers of events proceeding via l^Dj states (35.1) and via 
the 2^5*1 state (38.4) are approximately equal, with the dominant roles played by the 
transitions 3S ^ 2^ Pi ^ 1^D2 ^ l^Pi ^ IS (20.1 events) and 3S ^ 2^ Pi ^ 21^ Si ^ 
l^Pi IS (15.9 events). Equal numbers of 47/x+/i~ events are expected if muons can 
also be identified. The initial and final photon energies for these two sets of transitions 
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Table III: Predicted numbers of A'ye'^e events corresponding to 3S — > 2P ID 
IP ^ IS ^ e+e- or 3S ^ 2P ^ 2S ^ IP ^ IS ^ e+e' per 10^ T(35) decays. 
Numbers following spectroscopic symbols represent photon energies in MeV, in cm. of 
decaying state. 
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are the same, since both sets of transitions proceed via ^Pi states. However, the two 
intermediate pairs of energies are different: 99 and 261 MeV for the transitions via 
1^L>2, and 131 and 229 MeV for the transitions via 2 5*1. Thus, it should not be too 
hard to distinguish the two processes from one another. 

The branching ratio predictions of Ref. 0] were performed under the assumption 
that the hadronic widths of the D-wave states could be calculated purely using their 
color-singlet bb components. A more up-to-date calculation of hadronic widths, based 



on inclusion of color-octet bb contributions |^, is probably called for, but is beyond 



the scope of the present note. Even 100% augmentation of the hadronic widths of 
the l^Dj states would have little effect on their branching ratios for radiative decays, 
which are expected to be dominant. The partial widths for 2^Dj states to decay non- 
electromagnetically were not calculated in Ref. but were expected to be comparable 
to those for the l^Dj levels. 

We conclude with a discussion of energy scans in e+e^ annihilations for direct 
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production of bb{^Di) states. The results of Table indicate that these states are 
expected to lie between 6 and 11 MeV/c^ below the D-wave spin-weighted centers of 
gravity. Taking account of the predictions in Fig. |I|, one then expects the l^Di level 
to lie between 10.13 and 10.17 GeV/c^ while the 2^Di level should lie between 10.42 
and 10.46 GeV/c^. These predictions ignore coupled-channel distortions due to BB 
threshold [2^], as mentioned. 

As mentioned in Ref. 0, present limits on leptonic widths of the l^Di and 2^Di 
states p2| are about a factor of 10 to 15 above the predicted values of ree(l^-Di, 2^_Di) : 
(1.5, 2.7) eV. [Ref. [|ll|] finds Tee{l^Di) = 1.6 eV.] The CUSB Collaboration's search 
in the range from 10.34 to 10.52 GeV/c^ sets 90% c.l. upper limits of Tee{2^Di) < 40 
eV in this range on the basis of 5 pb~^ of integrated luminosity. Thus, an effective scan 
of this same range with 15 times the sensitivity would require at least 1 fb~^. Similar 



estimates have been obtained by the CLEO Collaboration p3 . 

Once a ^Di state has been found in an energy scan, does it have any distinctive final 
states? For the l^Di state, the dominant decay, with 60% branching ratio, was found 
H to be 7 + l^Po, with = 285 MeV. The I^Pq is expected to decay 97% of the time 
to hadrons. However, its hadronic decays do not appear to have the expected signature 
of a pair of nearly back-to-back jets p4|. For the 2^Di state, the corresponding photon 
energy for the dominant final state 7 + 2^Po is predicted to be 202 MeV. The 
branching ratio to 7 + l^Po, leading to a 559 MeV photon, is expected to be about a 
factor of 5 lower. 

We thank R. S. Galik and T. Skwarnicki for asking questions which led to this 
investigation, and for extensive discussions. We are grateful to C. T. H. Davies for 
communications with regard to lattice gauge theory results. This work was supported 
in part by the United States Department of Energy through Grant No. DE FG02 
90ER40560 and the Natural Sciences and Engineering Research Council of Canada. 
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